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1 INTRODUCTION 



ABSTRACT 

We present results and applications of high precision timing measurements of the 
binary millisecond pulsar J1012+5307. Combining our radio timing measurements 
with results based on optical observations, we derive complete 3-D velocity information 
for this system. Correcting for Doppler effects, we derive the intrinsic spin parameters 
of this pulsar and a characteristic age of 8.6 ± 1.9 Gyr. Our upper limit for the orbital 
eccentricity of only 8 x 10"'' (68% C.L.) is the smallest ever measured for a binary 
system. We demonstrate that this makes the pulsar an ideal laboratory to test certain 
aspects of alternative theories of gravitation. Our precision measurements suggest 
deviations from a simple pulsar spin-down timing model, which are consistent with 
timing noise and the extrapolation of the known behaviour of slowly rotating pulsars. 

Key words: binaries: general — pulsars: general — pulsars: individual: J1012-f 5307 
— gravitation — relativity — time 
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scope at Jodrell Bank. Nicastro et al. (1995) showed that the 
pulsar is in a binary system with an orbital period of 14.5 hr 
and a companion mass between 0.11 and 0.25 M© (90% 
confidence level, C.L.). Optical observations reported by 



Lorimer et al. (1995) revealed an optical counterpart within 
0'.'2±0'.'5 of the pulsar timing position, consistent with its 
being a white dwarf (WD) companion of mass 0.15 M0. 

Indeed, PSR J1012-I-5307 is one of the few examples of 
binary millisecond pulsars where optical observations of the 



mate s to correct for apparent acceleration effects ( Bhklovskii 
197C). A proper motion obtained by pulsar timing is usually 
converted into a transverse velocity by using a distance esti- 
mate derived from the dispersion in the interstellar medium, 
applying a model of th e free electron density distribution 
(Taylor & Cordes 1993). However, this model yields a high 
level of systematic uncertainty (e.g. Toscano et al. 1999b). 
The dispersion measure (DM) distance of 520 pc is indeed 
somewhat smaller than the distance derived from luminosity 
models. 



Applying classical cooling models (Iben & Tutukov 



WD CO npanion considerably enhance our knowledge about 



the radi o pulsar and its bi nary systemT^an Kerkwijk et 
al. 1996). Callanan et al. (1998), for instance, used their 



results from optical observations to obtain a number of es- 
sential pieces of information for this system. Comparing the 
measured optical luminosity to its value expected from WD 
models, they determined the distance to be d = 840±90 
pc. The Doppler shift of the measured H-spectrum of the 
WD companion gives an additional radial velocity compo- 
nent of 44±8 kms"^ relative to the solar system barycentre, 
SSB. From the radial velocity (semi-)amplitude and the or- 
bital parameters known from pulsar timing, the mass ratio 
between the pulsar and its companion is measured to be 
mpsR/rric =10.5±0.5. 

In order to derive the intrinsic spin-down age of mil- 



1986) for low mass WDs and the DM distance, Lorimer et 
al. derived an age of only 0.3 Gyr for the WD companion 
which was in apparent contradiction with the inferred spin- 
down age of the pulsar of 7 Gyr. However, Alberts et al. 
( [L996| ) were able to model the binary parameters if they as- 
sumed that the cooling age of the WD is similar to the spin- 
down age of the pulsar. More recently, Driebe et al. (199J ) 



modeled the cooling process of WDs, obtaining an age of 
6±1 Gyr for the companion of PSR J1012-I-5307. Ergma et 
al. (1998) modeled evolutionary sequences of neutron stars 



with close low mass binary comp anions and short orbital pe- 
riods, whilst Sarna et al. ( 1998 ) developed an evolutionary 
scenario for PSR J1012-I-5307, reproducing the observed bi- 
nary period and companion mass. They modeled the age of 
the neutron star to be between 4.5 and 6.0 Gyr, whilst the 
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mass of the WD was found to be 0.16 ± 0.02 M©, consis- 
tent with mass estimates derived by other mod els based on 
different observations and e volutionary models (Callanan et 
al. 199^; [Priebe et al. 1998[ ). 

In this paper, we combine a set of high precision tim- 
ing data with the information derived for the binary system 
from radio pulsar and optical WD observations as well as 
from model calculations. After a short description of the tim- 
ing observations and data reduction techniques, we describe 
a new binary timing model that is important for the low 
eccentricity case, discuss the resulting solution using statis- 
tical tests, and present the derived timing parameters. These 
results are analysed in the following sections. We show that 
proper motion and Galactic acceleration of PSR J1012-I-5307 
have a significant influence on the derived characteristic age. 
Furthermore, we obtain strong limits on the true eccentricity 
of the system and discuss its meaning for evolutionary sce- 
narios. As a result of the stringent limits obtained both on 
the change in orbital period as well as on the true eccentric- 
ity of the system, we demonstrate that this binary system is 
highly suitable for testing different theories of gravitation. 



2 OBSERVATIONS 
2.1 Effelsberg timing 

We have made regular high precision timing observations of 
PSR J1012-I-5307 since October 1996 using the 100-m radio 
telescope of the Max Planck Institut fiir Radioastronomie in 
Effelsberg near Bonn. The typical observing rate has been 
once per month with a usual observing time of about one 
hour. The overall post-fit RMS of the timing model ap- 
plied to the 1213 resulting times-of-arrival (TOAs) is 3.1 
^s. Sometimes, however, the source shows strong intensity 
fluctuations on time scales of typically two hours caused by 
scintillation effects in the ionised component of the inter- 
stellar medium. We made use of these intensity maxima by 
optimising our observing strategy, i.e., staying "on source" 
until the pulsar becomes weaker again. As a consequence of 
this highly successful "scintillation hopping", many TOAs 
are of much higher quality, with errors frequently smaller 
than 1 /is. Figure |l| shows a typical example for the develop- 
ment of TOA uncertainties during a scintillation maximum. 
Since we have many such high quality TOAs covering sev- 
eral orbits, they proved to be extremely useful for the precise 
determination of orbital parameters. 
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Figure 1. Errors of the TOAs during a scintillation maximum. 



the primary focus with a bandwidth of 30 MHz. The system 
temperature during the observations was about 65 K on cold 
sky with a telescope gain of about 1.5 K Jy~^. The 2700- 
MHz data were obtained with a cooled HEMT receiver in 
the secondary focus, with a system temperature of typically 
50 K on cold sky. As for the other receivers, LHC and RHC 
signals were mixed down to an IF of 150 MHz. 



2.1.2 Coherent De-Disperser 

The signals received from the telescope were acquired and 
processed with the Effelsberg-Berkeley Pulsar Processor 
(EBPP) which removes the dispersive effects of the inter- 
stellar medi um on-line using the techn ique of "coherent de- 
dispersion" (Hankins & Rickett 1975). Before entering the 



EBPP, the two IF LHC and RHC signals are converted to an 
internal IF of 440 MHz. A maximum bandpass of 112 MHzQ 
is split into four independent portions for both circular po- 
larisations. Each portion is mixed down to baseband and 
subdi vided into 8 narrow channels via a set of digital fll- 
ters (Backer et al. 1997) for coherent online de-dispersion 



by convolution using a custom VLSI device. In total 64 out- 
put signals are detected and integrated in phase with the 
predicted topocentric pulse period. 



2.1.1 Receiving systems 

Most of the data were obtained with a 1.3-1.7 GHz tun- 
able HEMT receiver installed in the primary focus of the 
telescope. The noise temperature of this system is 25 K, re- 
sulting in a system temperature from 30 to 50 K on cold sky 
depending on elevation. The antenna gain at these frequen- 
cies is 1.5 K Jy~^. An intermediate frequency (IF) centred 
on 150 MHz for left-hand (LHC) and right-hand (RHC) cir- 
cularly polarised signals was obtained after down-conversion 
from a central RF frequency of usually 1410 MHz. 

In order to monitor DM changes, occasionally we also 
obtained data at 860 and 2700 MHz. For the 860-MHz ob- 
servations we used an uncooled HEMT-receiver operating in 



2.1.3 Arrival times measurements 

A TOA was calculated for each average profile obtained dur- 
ing a 5-10 min observation. During this process, the ob- 
served time-stamped profile was compared to a synthetic 
template, which was constructed out of 12 Gaussian com- 
ponents fitted to a high signal-to-noise standard profile (see 
Kramer et al. 1998; 1999). This template matching was done 
by a least-sq uares fitting of the Fourier-transformed data 
( [Taylor 19911 ). The final TOA was obtained by using the 



* The bandwidth available depends on observing frequency and 
DM of the pulsar. At 1.41 GHz we used a bandwidth of 56 MHz 
for PSR J1012-I-5307. 
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measured time delay between the actual profile and the tem- 
plate, and the accurate time stamp of the data provided by 
a local H-maser and corrected off-line to UTC(NIST) using 
recorded information from the satellites of the Global Po- 
sitioning System (GPS). The uncertainty of each TOA was 
estima ted using a method described by Downs & Reichley 
(I1983D. 



2.2 Jodrell Bank timing 

Since its discovery in 1993, PSR J1012-I-5307 has been reg- 
ularly monitored using the 76-m Lovell telescope at Jodrell 
Bank with cryogenic receivers at 408, 606, and 1404 MHz. 
Both LHC and RHC signals were observed using a 2 x 64 x 
0.125-MHz fiher bank at 408 and 606 MHz and a 2 x 32 x 1.0- 
MHz filter bank at 1404 MHz. After detection, the signals 
from the two polarisations were filtered, digitised at appro- 
priate sampling intervals, and incoherently de-dispersed in 
hardware before being folded on-line with the topocentric 
pulse period and written to disk. Each integration was typ- 
ically of 1-3 minutes duration; 6 or 12 such integrations 
usually constituted an observation. In the analysis stage, 
the profiles were added in polarisation pairs before being 
summed to produce a single total-intensity profile. A stan- 
dard pulse template was fitted to the observed profiles at 
each frequency to determine a total of 400 pulse TOAs. De- 
tails of the observing system and the data reduction scheme 
can be found elsewhere (e.g. Bell et al. 1997). 

Although the templates used in the Effelsberg and Jo- 
drell Bank data reduction differed, resulting offsets were ab- 
sorbed in a global least-squares fit. The typical measurement 
accuracy of the Jodrell Bank 1404 MHz timing data is about 
9 /is. 



3 PRECISION TIMING 

The combined TOAs, corrected to UTC(NIST) and 
weighted by their individual uncertainties determined in the 
fitting i jirocess, were analysed with the tempo software pack- 
ag^, usi ng the DE200 e phemeris of the Jet Propulsion Lab- 
oratory (Standish 1990). Owing to the extremely low eccen- 
tricity of this system, we applied a binary timing model that 
uses the Laplace Lagrange parameters r] and k and the time 
of ascending node Tasc in place of the Keplerian parame- 
ters e, ijj and To (see Appendix A for details). The Keplerian 
parameters as a function of those used by the model can be 
calculated as 



e 
To 



arctan(r;/K) 

Tasc + ^ arctan(r7/ k) 
In 



(1) 
(2) 

(3) 



TEMPO minimizes the sum of weighted squared timing 
residuals, i.e., the difference between observed and model 
TOAs, ri, yielding a set of improved pulsar parameters 
and post-fit timing residuals. The timing residuals shown 
in Fig. ^ indicate deviations from the timing model which 
are slightly larger than the calculated TOA uncertainties, 

t http:/ /pulsar. princeton.edu/tempo 



(Ti. This most likely stems from a small underestimation of 
the error in the TOAs. The data uncertainties used in the 
following were scaled by an appropriate factor to achieve a 
uniform reduced — 1 for each data set. In order to test 
the statistical properties of our data we have generated his- 
tograms of the deviations of the timing residuals from the 
model. The final residuals show the expected Gaussian dis- 
tribution (Fig. 3). 

In order to search for systematic errors that are either 
small enough to be individually undetectable or have a time 
dependence which approximates some linear combination of 
terms in the model, we have searched for tell-tale correla- 
tions among nearly adjacent post-fit residuals (cf. Taylor & 
Weisberg 1989). We combined 2, 4, 8 etc. consecutive TOAs 
and evaluated the average residual of the remaining data 
set. A "clean" data set should exhibit a slope of —0.5 in a 
plot of RMS versus number of combined TOAs. The results, 
as shown in Fig. ^, show no or very little correlation among 
Effelsberg data and the Jodrell Bank data recorded at 410 
MHz after employing the timing model described in Table ^ 
The JodreU Bank data at 1400 MHz and 606 MHz exhibit, 
however, some correlation. We tested for possible impact of 
these correlations on the fit parameters by analysing the 
timing data without one or both of these data sub-sets. The 
parameters derived from these reduced data sets did not 
significantly differ from those derived by the fit to all the 
data. 

The final parameters obtained from fitting the timing 
model to our data of PSR J1012-I-5307 are shown in Ta- 
ble |l]. The upper limits on the derivatives of e, x and Pb 
were obtained by separately fitting for each of them. These 
were set to zero during the fits for all other parameters. The 
residuals show slow variations that can be removed by fit- 
ting for P. We do not believe that instrumental effects can 
account for this effect but rather suggest that the measure- 
ment reflects rotational instabilities in the pulsar. In this 
case, PSR J1012-I-5307 is only the second Galactic field mil- 
lisecond pulsar to show a significant second spi n frequency 



derivative, i>, the other being PSR B1937-I-21 (Kaspi et al 



1994| )p'|. We will discuss the possible important imphcations 
of this result in Section ^ 



4 INTRINSIC PULSAR PARAMETERS 

As summarised in Section ^ the binary millisecond pulsar 
PSR J1012-I-5307 is unique since the combination of optical 
and radio measurements allows a full determination of the 
motion of the system through space and also the determi- 
nation of several intrinsic parameters. 



4.1 3-D velocity 

The full 3-d motion of the pulsar relative to the SSB is 
determined by combining the timing proper motion with 
distance to the system obtained from luminosity models for 



"f Cognard et al. (1996) measured a u for the millisecond pulsar 
B1821— 24 located in a globular cluster. Here the acceleration 
by the cluster gravitational potential may significantly affect the 
observed spin-down parameters. 
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Right Ascension, a {J2000) 
Declination, 6 (J2000) 
^ic (mas yr-l) 
^lS (mas yr^i) 

V (Hz) 
r> (s-2) 
i> (s-3) 

P (ms) 
P 

P (s-l) 
Epoch (MJD) 



10'U2™33!43370(3) 

53°07'02'.'5884(4) 

2.4(2) 

-25.2(2) 

190.267837621903(4) 
-6.2029(4) ■ 10-1*^ 
-9.8(2.1) • 10-27 

5.2557490141197(1) 
1.7134(1) ■ 10-2" 
5.1(1.1) ■ 10-3^ 

50700.0 





Dispersion Measure (DM) (cm-^' pc ) 9.0233(2) 



Figure 3. Distribution of all post-fit timing residuals normalized 
by their uncertainties. 



Orbital period, Pj, (days) 
Proj. semi-major axis, x (s) 

■n 

Tasc (MJD) 



Upper limits: 

Parallax, 7r(mas) 
e 

e (s-l) 



DM (cm 3 pc yr 1) 



0.60467271355(8) 
0.5818172(2) 
0.9(8) ■ 10-"^ 
0.01(80) ■ 10-8 
50700.0816290(1) 



< 1.3 

< 1.3 X lO-'' " 

< 2 ■ 10-1* 

< 1.4 ■ 10-1* 

< 1 ■ 10-13 

< 1.2 ■ 10-* 



° See section [4.4| for details. 

Table 1. Timing parameters for the millisecond pulsar 
J1012-I-5307. Quoted errors correspond to twice the errors de- 
rived by TEMPO. Upper limits represent 95% C.L. 
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Figure 2. Post-fit timing residuals as a function of observing 
year. 



Nb. of consecutive residuals averaged 

Figure 4. Mean post-fit RMS residuals of the data sets versus 
the number of consecutive averaged timing residuals. The symbols 
represent data from: EBPP 1.4 GHz (A) and 2.7 GHz (□), and 
Jodrell Bank 410 MHz (Q), 606 MHz (0), and 1.4 GHz (*). The 
870 MHz data from Effelsberg are not included here owing to 
their small number. 



d = 



-102 ± 11 kms" 



(5) 



Optical measurements (cf. Section |l|) yield a radial velocity 
component of -1-44(8) kms~i. Hence, the total velocity of the 
system is 111(28) kms-i. This value is consistent with the 
average space velocity of millis econd pulsars of 130 km s- i 
estimated by Lyne et al. (1998) and Toscano et al. (1999a). 



the white dwarf companion, d — 840 ± 90 pc. This distance 
is consistent with the obtained upper limit for the timing 
parallax. Given the pulsar's ecliptic latitude of 38°, it is not 
unexpected that a significant value for the timing parallax 
could not be obtained. We derive transverse velocities of 



• d 



8.8 ± 1.2 kms" 



(4) 



4.2 Doppler effects 

Due to the motion of PSR J1012-I-5307 relative to the solar 
system, any period P°^^ of the pulsar spin or binary sys- 
tem, measured at the SSB, differs from the period P'"* in 
the pulsar reference frame by a Doppler factor D. FoUow- 
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ing Damour & Taylor (1991), we write the relation between 
these periods 



pobs _ pintr _ ^ 



(6) 



with n the unit vector to the pulsar a nd Vr the relative 



velo city between the pulsar and the SSB ( Damour fc Taylor 



1992 |). | \s Vr/c is typically less than 0.1% for millisecond 
pulsars, this Doppler shift is for most purposes unimportant. 
Hence, we will for the rest of this paper denote P°^^ as 
P. However, the period derivative may be modified by the 
relative motion. 

Calculating the time-derivative of Eqn. ^and separating 
the effe ct of the proper motion of the system (Shklovskii 
om the influence of the Galactic acceleration, we 
derive as contributions to the period derivative 



pintr jj _^ pShk ^ pG 



(7) 



where the difference of the constant Doppler factor D from 
unity may be neglected. The contributions P^*""^ + P^''' = 
P° are written explicitly 



i)Shk 1 2 , o 

— — a dF 

c 



(8) 
(9) 



where d is the distance to the pulsar. Here npsr is the posi- 
tion vector to the pulsar, fi its proper motion and apsr and 
ao the Galactic acceleration of the pulsar and the Earth, 
respectively. We measure npsr and /i from our timing obser- 
vations. The small correction term p'^^^/p = — 7 x 10^^"s~^ 
for the Galactic acceleratio n is obtained applving a rnodol 



Innanen 1987 



Kui- 



for the Galactic pote ntial ([Carlberg^ 

ken & Qilmore 1989). Inserting these parameters into Eqns 



and ^ , we obtain a Doppler correction of 
P°/P = 1.4 ±0.3x10" 



(10) 



Its error is dominated by the uncertainty in the distance to 
the pulsar. 



4.3 Characteristic age 

We can now use our timing results to obtain the 
true spin-down age, the so-called characteristic age of 
PSR J1012-H5307. With a spin period of P = 5.256 ms and 
a Doppler factor as given in Eqn. |lo| we derive a Doppler 
correction of 



7.4 ± 1.6 X 10" 



(11) 



which has to be subtracted from the measured P before 
calculating the intrinsic characteristic age of the system: 



2 pint 2(P°''=-PD) 



8.6± 1.9Gyr. 



(12) 



We note that the characteristic age assumes magnetic dipole 
braking (i.e., braking index of 3) and is only a realistic es- 
timate of the real pulsar age if the initial spin period was 
sufficiently small (Camilo et al. 1994). Although this may 
not apply to He- WD binaries (Backer 1998), the character- 
istic age is in reasonable agreement with a coolin g age of the 
white dwarf as derived by evolutionarv mode ls ( Alberts et 
al. 1996|; |Sarna et al. 1998|; [Priebe et al. 199^). 




1 1.5 

eccentricity [10 

Figure 5. Distribution of values for the real eccentricity from 
simulated parameters (see text). 



4.4 Orbital eccentricity 

In Appendix A we analyse the contribution of the Shapiro 
delay for binary orbits with moderate inclination angles and 
negligible intrinsic eccentricities to the apparent eccentric- 
ity of the system. As we show, the Shapiro delay cannot be 
separated from the Roemer delay, which leads to a small 
correction to the observed eccentricity. In order to investi- 
gate this effect for PSR J1012-I-5307 quantitatively and to 
obtain the true orbital eccentricity, we make use of indepen- 
dent determinations of the system parameters. 

From the companion mass of 0.16(2) Mq, the measured 
mass ratio of g = mp/rric = 10.5(5), and the mass function 



{nic + nip) 



- 0.000578 M, 







(13) 



obtained from the timing analysis, we derive the range r and 
shape s of the Shapiro delay in the system according to 



sm t = 



4.9255 (mc/Mo) and 



(14) 
(15) 



The resulting shape parameter of s « 0.8 indicates a moder- 
ate inclination of the system, i.e., i = 52° or 128°. The con- 
tribution of the Sh apiro delay to the observed rj is calculated 
according to Eqn. A22 and subtracted from the measured 
value. Applying Eqn. ^, we obtain upper limits on the true 
eccentricity of 



e'"*' < 0.8 X 10" 



e'"*" < 1.3 X 10" 



70 C.L.) 

(95% C.L.), 



(16) 



which were derived from Monte Carlo simulations. In these 
calculations, whose results are displayed in Fig. ^, we simu- 
lated the parameters q and rric and the observed k and ry in 
accordance with the observational uncertainties. 



5 TESTS OF THEORIES OF GRAVITATION 
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5.1 Dipole gravitational waves 

Unlike general relativity, many alternative theories of grav- 
ity predict the presence of a radiative monopole and d ipole 



as well a s the quadrupole and higher order multipoles (Will 
1993 |). ^ calar-tensor theories, for instance, predict for bi- 
nary systems a loss of orbital energy which is at highest 
order dominated by scalar dipole radiation. As a result, the 
period Pb of a circular orbit will change with 



dipole 



47r G, nip rUc 
c3 Pi, nip + nic 



(Qp ■ 



(17) 



where G, is the bare gravitational constant, c is the speed 
of light and m„ and rric the mass of the p ulsar and its com- 
panion ( Damour fc Esposito-Farese 199q ). The parameters 
Up and Qc represent the effective coupling strengths of the 
scalar field to the pulsar and its comp anion, respectively. 

Damour & Esposito-Farese (1996) show that can be 
neglected compared to Qp, since the gravitational binding 
energy per unit mass of a WD, £wd = E^''^^ /mc? « 10"'*, 
is very much smaller than that of a 1.4-Mq neutron star, 
Ep ^ 0.15. 

Using this assumption, and taking G* ~ G, Eqn. ^ can 
be approximated as: 



pdipolc 



nic 



q 2 



(18) 



Since all terms except Op in this equation are measured, 
the effective coupling strength of PSR J1012+5307 can be 
restricted by limits on Pb- We will show that these are the 
tightest bounds on Qp ever measured for a neutron star. 

The observable rate of change of the orbital period of 
PSR J1012-I-5307 is the sum of several contributions, either 
intrinsic to the orbit or caused by projection effects: 



pG _|_ pdipolc 



(19) 



(Damour & Taylor 1991). Here, the Doppler correction of Pb 
is given by and the contribution -Pj"^ represents gravita- 
tional radiation as predicted by general relativity. The values 
of these effects can be estimated and may be subtracted from 
the measured value of Pb- The additional contributions P^ 
and p^'*'P°'° represent the period derivative due to a change 
of the gravitation constant and due to gravitational dipole 
radiation, respectively. All terms of gravitational radiation 
not accounted for are of the order or higher. We derive 
PfP from Eqn. |l^, with the orbital period of 52 240 s, to be 

P(P = 7.3± 1.6 X 10"^*. (20) 

The orbital period decay derived from the theory of general 
relativity for a circular orbit can be approximated to 



384 TT^ (Gmcf/^ 



pV3^5 (g 4-1)1/3 



-1.0±0.2xl0' 



,(21) 



using the WD-mass rric = 0.16(2) Mq and t he mass ratio of 
q ~ 10. 5 between the pulsar and the WD ( Callanan et al 
199f~| 

A variation in the orbital period as a result of a possible 
time dependence of the gravitational constant can for this 
system be approximated to 



-2|(l + Cp)P, 



(22) 



(Nordtvedt 1990) where Cp is the compactness of the pulsar, 
which is approximately two times its b inding energy. Using 
the upper limit \G/G\ <, 8x 10"^Vr" 
we can restrict 



^ (Williams etal. 1996) 



jPfci < 3.6 X 10-^ 



(23) 



This value is only about 1/3 of our 1-a error for Pb. We 
can therefore exclude a significant influence on Pb from this 
effect. 

Subtracting the known contributions from the limit 
measured for Pb yields a value of p^'P°i'= ~ —0.6(1.1) 10"^^. 
Using this value we obtain limits for ap of 



Qp < 2 X lO""* 
ftp < 4 X 10"* 



(68% C.L.) 
(95% C.L.), 



(24) 



where we have taken into account uncertainties in the pa- 
rameter values used. 

Simulations of future observations show that the preci- 
sion of the measurement of Pi, will increase significantly. 
A continuation of our timing observations with an accu- 
racy and observing rate similar to that of the past few 
years will lead to an improvement by a factor of about 3 
within the next three years. The applicability of this im- 
provement for the value of ap might be reduced by the un- 
known contribution of G, unless the limits on this parameter 
are also improved. It is also desirable that the distance to 
PSR J1012-I-5307 be determined with higher accuracy to al- 
low a more reliable estimate of the transverse Doppler effect. 

5.2 Local Lorentz invariance 

Will & Nordtvedt (1972) pointed out that one expects the 
existence of a preferred cosmic rest frame for gravitational 
interaction if gravity is mediated in part by a long range vec- 
tor field or by a second tensor field. In the post-Newtonian 
limit all the gravitational effects associated with such a pre- 
ferred frame are phenomenologically describable by two pa- 
rameters, ai and Q2. Only in gravity theories without a pre- 
ferred frame, e.g. general relativity, are the values of these 
parameters equal to zero. 

A consequence of a preferred cosmic rest frame would 
be the introduction of an eccentricity into orbits of gravita- 
tionally bound bodies of different masses which are in mo- 
tion relative to this preferred fram e. This effect is quantified 
by the parameter qi. Wex (200C) has applied a statistical 
analysis of all pulsar binary systems to find limits for this pa- 
rameter. This analysis yields |ai| < 1.4 x lO""* at 95% C.L. 
Due to its extremely low true eccentricity and the known 
three dimensional proper motion, the data presented here 
for PSR J1012-f 5307 and included in Wex's analysis play a 
central role in this analysis. They might provide more strin- 
gent tests of local Lorentz invariance, if the real eccentricity 
of the system is even smaller and better observational limits 
on it are found. 



6 DISCUSSION 
6.1 Timing noise 

Over the observed time span of about 7 years, 
PSR J1012-I-5307 shows a significant deviation from a sim- 
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pie v-O spin down behaviour. In order to model this effect, 
the spin-down model can be extended by allowing for a non- 
zero second derivative of the spin frequency, i> (e.g. Cordes 
& Downes 1985). However, rather than being the result of 
rotational irregularities of the neutron star itself, a possible 
alternative explanation for a second spin frequency deriva- 
tive could be a third order term of the proper motion. 



..D o 2 

V = — 3 I'intr 

C 



2 0- 

^ I'mtr 

C 



(25) 



where Vr is the radial velocity of the system (Phinney 1992). 
This term can be computed and it is of order 10"'^" s~'^, i.e., 
too small to explain the measured effect. 

In order to investigate the possibility that the observed 
timing noise is caused by small changes in the DM, we mod- 
eled our multi-frequency data with a simultaneous fit of v 
and a second-order polynomial for DM. While DM and DM 
were not significant, the u measurement was still a 5-(t ef- 
fect. The upper limit on DM given in Table ^ is interestingly 
much smaller than one would expect for DM=9.02 cm~^ pc 
(Backer et al. 1993). 

As a possible explanation, a second period derivative in 
the timing residuals can also be the result of a sec ond light 
companion in a wider orbit (Thorsett et al. 199!:). Assum- 
ing a circular orbit and an intermediate orbital phase, the 
measured value for D can be explained with e.g. a planet of 
terrestrial mass at ~ 30 AU distance or a Jupiter-like planet 
at ~ 170 AU distance from the centre of mass. Future ob- 
servations will show if the timing residuals can be explained 
by a second Keplerian orbit. 

However, it is more likely that the irregularities are in- 
trinsic to the pulsar itself. Timing noise has been stud ied fo r 
normal, slowly rotating pulsars. Arzoumanian et al. (1994), 
for instance, defined the stability parameter 



A(t) = log(^t3 

D V 



(26) 



where t is the observing span which in their case was 10 
s. Despite noticeable scatter, they obtained a correlation of 
this parameter with the period derivative for slowly rotating 
pulsars. This correlation predicts a value of Ag ~ —5.3 for 
PSR J1012-I-5307. Using timing data over a time span 10* s, 
we measure a value of Ag ~ — 5. Given that the fit leading to 
Arzoumanian et al.'s relation was made by eye, both values 
are in excellent agreement, in particular given the large scat- 
ter seen for normal pulsars. Moreover, the value of the sec- 
ond spin frequency derivative, v — —9.8(2.1) x 10~^^ s~^, is 



similar to t he value of v published for PSR B1937-I-21 (Kaspi 
et al. 1)94). Both sources exhibit a Ag consistent with the 
correlation derived by Arzoumanian et al. Although we can- 
not yet completely rule out a different origin of the measured 
period second derivative, the value determined is apparently 
consistent with what is known for timing noise from normal 
pulsars. Assuming that we are now indeed observing rota- 
tional irregularities in two Galactic field millisecond pulsars, 
being consistent with the behaviour of normal pulsars ex- 
trapolated to very small period derivatives, it may indicate 
that the empirical relation found by Arzoumanian et al. is 
generally valid for millisecond pulsars as well as for normal 
pulsars. As other millisecond pulsars are now being moni- 
tored with similar timing precision over similar time spans, 
this conjecture can be tested. It may well be the case that 




000 



iinary Period (days) 



Figure 6. Relation between binary period and orbital eccentric- 
ity (filled symbols are measurements, open symbols upper limits) 
in binary systems with a millisecond pulsar and white dwarf com- 
panion. 



we are finally exploring the ultimate accuracy of the pul- 
sar clockwork, deciding over the suitability of millisecond 
pulsars as potential standards of time. 



6.2 Orbital eccentricity 
evolution 



a relic of binary 



The c onvective fluctuation-dissipation theory of Phinney 
( [l992| ) predicts a strong correlation between the orbital ec- 
centricities and orbital periods of binary pulsars that have 
been recycled by stable mass-transfer from a Roche-lobe- 
filling, low-mass red giant. The measurements and upper 
limits for the orbital eccentricities of systems with appro- 
priate evolutionary history are plotted in Fig.|^. They are, 
with the exception of PSR J0613— 0200, in excellent agree- 
ment with the model of Phinney & Kulkarni (1994; curves 
in Fig. |). 

The lowest orbital eccentricity of a millisecond binary 
pulsar system known is the up per li mit of 1.3xl0~^ (68% 
C.L.) derived by Camilo et al. ( |l996| ) for PSR J2317+1439, 
which has an orbital period of 2.5 days. While this source 
follows the predicted trend, the eccentricity of the 1.2- 
d binary pulsar J0613— 0200 has been determined to be 
3.8 ± 1.0 X 10"^, which exceeds the 95% C.L. of this model 
relation by a fac tor o f about 2.5 (Fig. ^). However, Phin- 
ney & Kulkarni (1994) pointed out that for extremely close 



binary systems the orbital eccentricity is difficult to model. 
Therefore, it is interesting to note that the upper limit for 
the eccentricity of the 0.6-d binary pulsar J1012-I-5307 of 
only e < 0.8 X 10"^ (68% C.L.) is significantly lower than 
that of PSR J0613-0200. 



6.3 Future observations 

The motion of the binary system relative to the SSB results 
in a change of the orbital inclination that causes a variation 
in the projected semi-major ax i s of the system (lArzouma- 



nian et al. 1996; Kopeikin 1996; BeU et al. 1997 



andhu et 



al. 1997). The amplitude of this effect is given by 
= 1.54 X lO'^'^x cot i {—jia sin f2 + cos Q) 



(27) 
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where SI is the position angle of the ascending node and 
and the components of the proper motion in masyr"^. 
Given the value of the orbital inclination sinz ~ 0.8 as de- 
termined from the white dwarf optical light curve, a mea- 
surement of X allows one to restrict the orbital orientation 
Q.. Using our current data set, we derive a 1-a uncertainty 
for X that is about twice the value expected for a perfect 
alignment between and ^. Continuing timing observation 
will permit the measurement of this effect. Simulations, as- 
suming an accuracy and data-taking rate similar to those at 
present, show that a 3-cr detection of x should be possible 
within the next four or five years. This would finally com- 
plete our knowledge of the full orientation and motion of the 
PSR J1012-(-5307 system. 



7 SUMMARY 

We have performed long-term high-precision timing obser- 
vations of PSR J1012-f 5307. Combining our measurements 
with results based on optical observations, we derived com- 
plete 3-D velocity information for this system, permitting 
the correction of the measured spin parameters for Doppler 
effects. Due to the precision of our measurements and the 
extremely small eccentricity of this binary pulsar, we could 
use it as an ideal laboratory to test certain aspects of the- 
ories of gravitation. The deviation of the timing data from 
a simple spin-down model is probably due to rotational in- 
stabilities which are consistent with the extrapolation of the 
known spin-down behaviour of slowly rotating pulsars, sug- 
gesting that this phenomenon is a common property of most 
pulsars. 
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APPENDIX A: TIMING OF 
SMALL-ECCENTRICITY BINARY PULSARS 

Neglecting relativistic effects and parallax, the barycentric 
arrival time, ft, of a signal emitted by a pulsar in an orbit 
around a compact companion is given by 



tft-to = r + Afl(r). 



(Al) 



T denotes the time of emission of the pulse signal as mea- 
sured in the reference frame of the pulsar, i.e., it is directly 
related to the rotational phase (j) of the pulsar by 

ct) = ct)Q + vT+'^vT'^ + ... (A2) 

The Roemer delay, A^, caused by the orbital m otion of the 
pulsar is given by ( iBlandford fc Teukolskyl97e| ) 



Afl = a; [(cos (7 — e) sinoj + sin ?7(1 — e^)'^''^ cosoj] . (A3) 

where the eccentric anomaly, U , is determined by Kepler's 
equation 



U ~ esmU = m (T~To) 



rib 



27r/P6, 



(A4) 



and uj denotes the longitude of periastron. 

For small-eccentricity binary pulsars, however, the lo- 
cation of periastron is not a prominent feature in the TOAs. 
Therefore, modeling the timing data of small-eccentricity 
binary pulsars with equation ( AS ) leads to very high corre- 
lations between the parameters ui and Tq. As a result, uj and 
To show unacceptably large uncertainties in the estima- 
tion of these parameters. 

Al A timing model for small-eccentricity binary 
pulsars 

Neglecting terms of order the orbital motion X = 
i?(cos (p, sin </5, 0) of a pulsar in a small-eccentricity binary 
system, is given by (see, e.g., Roy 1988) 



R = flp (1 - ecosM) 
ifi = M + 2e sin M 



M = nb{T-To), (A5) 



where M = (f = corresponds to the location of periastron. 
To first order approximation in the small eccentricity e the 
Roemer delay (Eq. A3) can therefore be written as 



(A6) 



Ar ~ X ^sin $ + ^ sin 2$ — ^ cos 2$ 

$ = 716 (r-r,3c), 

where terms which are constant in time are omitted. The 
three Keplerian parameters To, e, and u are replaced by the 
time of ascending node which is defined by 



T 

J- PIS 



To - cj/ub, 



(A7) 



and the first and second Laplace- Lagrange parameter, 

r] = e sinoj and k = e costj, (A-8) 

respectively. Note that the actual time when the pulsar 
passes through the ascending node is given by Tasc + 2ri/nb. 
The time of conjunction is given by Tasc + Pb/4: — 2k, /rib. 

Eq. (AC) accounts only for first-order corrections i n e . 
Therefor e, th e difference between the exact expression (A3) 
and Eq. ( A6) can grow up to xe^ . For most of the low eccen- 
tricity binary pulsars, the error in the TOA measurements 



is much larger than xe^ and thus the linear-in-e model is 
sufficient. 

Secular changes in the parameters rib, x, e, and u in 
Eq. (A3) are accounted for in the new timing model by using 

Ub (T — Tasc) + \nb [T — Tasc), 



and 



X = a;o-f i (T-Tasc), 

r) = r]0 +V {T — Tasc), 
K 



KQ + k {T — Tasc), 



in Eq. (A6) where the relations 



Ub 
Tasc 

hb 

V 
k 

hold. 



Ub +Lb — nb{To — Tasc), 



To 



rib, 



nb + UJ 



e smuj + e cos lo lo, 
e cos ui — e sin u to 



(A9) 



(AlO) 



(All) 
(A12) 

(A13) 
(A14) 
(A15) 



A2 Small-eccentricity orbits and Shapiro delay 

For small-eccentricity binary pulsars, the Shapiro delay can 
be written as 



-2r ln(l — ssin$). 



(A16) 



wher e r = Gmcjc^ and s = sini. As a Fourier series 
Eq. (A16) takes the form 



As = 2r{ao -I- bi sin <I> — a2 cos 2$ -f 
where 



ao 



&i = 2 



0.2 = 



In 



1 + VT 



1 - V]^s2 



41 



(A17) 

(A18) 

(A19) 
(A20) 



<^ 1 (ne arly edge-on) are 
. . in Eq. (A17), significant. 



Only for orbits where \J\ — s^ 
higher harmonics, indicated as 
Otherwise, the Shapiro delay cannot be separated from the 
Roemer delay. Consequently the observed values for x and 
rj differ from their intrinsic values by 

2.{o''s) _ ^ _|_ 2rbi , 

j^iobs) _ r) + 4:ra2/x. 



as can be seen by a comparison of Eq. ( A17 ) with Eq. 
In terms of the Blandford-Teukolsky model, Eq. ( |A22[ ) 
means that the observed value of the eccentricity e and the 
observed value of the longitude of periastron lo are different 
from the intrinsic values of these parameters. 

This new model has been implemented in the tempo 
timing software as binary model ELLl. 
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